Alpha-alpha interaction with chiral two pion exchange and ^Be lifetime. 
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Assuming that alpha particles are described by a scalar-isoscalar chiral invariant heavy field it is 
shown how chiral symmetry determines the alpha-alpha interaction at long distances unambiguously 
from dimensional power counting of an effective chiral Lagrangean. The leading strong contribution 
is given by a two pion exchange potential which turns out to be attractive and singular at the origin, 
hence demanding renormalization. When *Be is treated as a resonance state a model independent 
correlation between the Q-factor and lifetime l/F for the decay into two alpha particles arises. 
For parameters compatible with potential model analyses of low energy na scattering it is found a 
Breit-Wigner width F — 4.3(3)eV very close to the experimental value, Fexp. = 5.57(25)eV. 
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Low energy nuclear reactions in the energy range of as- 
trophysical interest are generally extremely hard to mea- 
sure experimentally in the laboratory This applies 
in particular to aa scattering where ^Be is produced 
as a narrow resonance In this paper we approach 
this problem from the theoretical side in the context of 
the chiral symmetry]^, and Effective Field Theories 
(EFT) 0. Specifically, it will be shown how the life- 
time of *Be in alpha-decay may be accounted for in a 
model independent fashion by assuming an elementary 
field for the a-particle, and exploiting the spontaneous 
breakdown of chiral symmetry of QCD as well as the fact 
that the ^Be ground state lies right above the aa thresh- 
old. In addition to the standard Coulomb interaction, 
the scalar-isoscalar character of the a particle implies 
that the long distance strong interaction is dominated 
by two pion exchange (TPE) regardless on any specific 
internal structure. 

Pion exchange interactions between a particles have 
been treated in the past in a variety of ways. A resonat- 
ing group method approach was used in Ref. [i*] with 
an approximation for the TPE in the mid-range. For- 
ward dispersion relations for aa scattering have been dis- 
cussed [A- A folding model from a NN potential was used 
in Ref. [8| and the /-wave phase shift was computed in 
first order perturbation theory. An EFT description of 
narrow resonances has been undertaken for pure contact 
theories 9], i.e. without pions or Coulomb forces. 

The idea of associating elementary fields to nuclei at 
low energies is rather natural [lol |. The a particle is a 
^He nucleus with (J^,T) = (0+,0), charge Z^e = +2e 
and mass = 2{Mp -f M„) - B = 3727.37MeV, with 
binding energy B = 28.2957MeV. Because of the ten- 
sor force, the wave function for the ground state can be a 
positive-parity mixture of three ^5*0, six "^Pqi and five ^Dq 
orthogonal states [lH , the symmetric S-wave component 
being the dominant part of the wave function, with signif- 



icant D-wave and almost negligible P-wave contributions. 
Thus, as a starting point we associate a scalar-isoscalar 
Klein-Gordon charged field, a{x) to the ^He nucleus. Un- 
der charge conjugation the a-particle should transform 
into an anti-a-particle a{x) a{x) meaning that the 
field is non-hermitean. Further, under SU{2)ii® SU{2)l 
chiral transformations we assume a{x) to be invariant; 
candidates for chiral partners with the same _B = 4 
baryon number would be *H and *Li having both spin-2, 
{J^,T) = (2^,1), thus belonging to different Poincare 
group representations. The effective Lagrangean will in- 
clude pions ^] and a particles which being much heavier. 
Ma 3> TOtt J are better treated by transforming the Klein- 
Gordon field as a{x) — e~''^''°''"'^ay{x) with ay{x) the 
heavy field and v'^ a four- vector fulfilling v'^ — 1, elimi- 
nating the heavy mass term 12 , . Keeping the leading 
Ma term, the effective Lagrangean reads 
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where the pion field in the non-linear representation is 
written as a S'?7(2)-matrix, U = e^'^''^/^ , with r the 
isospin Pauli matrices, / the pion weak decay constant 
in the chiral limit / = 88MeV, x = ?"^/2 and (,) 
means trace in isospin space. Here g^, gi, 52 and A 
are dimensionless coupling constants which are not fixed 
by chiral symmetry. This Lagrangean is the analog of 
the Weinberg- Tomozawa Lagrangean and EFT exten- 
sions for ttN interactions JL3,] to the case of the ira 
system. Photons are included by standard minimal cou- 
pling 9^Q! — + D^^a — d^a + ZaeiA^a. For definiteness 
we take. Ma = 3727.3MeV, f ^ 92.4MeV and 

TO ^ = 138MeV. 

To estimate the couplings go, 01 and g2 we look first 
at low energy ira scattering (see [J] for a review) . From 
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FIG. 1: Diagrams contributing to the a — a potential at long 
distances: One-Photon Exchange (left) and Two-Pion Ex- 
change loop (right). Full lines represent the a. particle field, 
the dashed lines are pion fields, wiggly lines are photon fields. 
The full blobs are the 00:7 and the aa-Kix vertices respectively. 



Eq. H]) we get the 7r°(p7r) -f a — > 7r''(p^) + a. invariant 
amplitude at lowest tree level order 



J 77 



{90A ■ -Ptt - giml + g2V ■ p^v ■ p'^) ,(2) 



where a and b are the final and initial pion isospin states 
respectively. In the limit ^ m.^ LAB and CM 
coincide, so taking v'^ — (1,0) and the kinematics as 
Ptt = and p'^ = {cj,p') with uj = \/p^ + rn'^ the 

scattering amplitude is given by {J- ~ Stty^/) 



Ua{P-,^) = 



go{ml+p'^{l ~ cosd)) -ffim^ -Fffa^^ 
Ao + Bo/ + SAi/ cos 6* . . . (3) 



where in the second line the threshold parameters Aq, 
Bq and Ai have been introduced. From mesic tt~ He 
atoms one has = (-0.138 -I- i0.045)fm 0, [li| while 
Bo = -O.lSfm^ and Ax = (0.42 -I- i0.06)fm^ from for- 
ward dispersion relations [151] . Our description is not 
realistic concerning the direct comparison with data; 
Coulomb distorsion has been disregarded and treating 
the OL particle as elementary precludes pion absorption 
since real 50,1,2 imply ImA/ — 0. However, realistic 
calculations allow to switch off these effects [16| yield- 
ing the pure potential values, = — 0.091 (17)fm and 
AP°* = 1.058(144)fm^ We take Bl°^ = -0.2(2)fm3 
from [1^, the error being an educated guess, since there 
are unresolved discrepancies (see e.g. [l7|). The small 
scattering length supports our perturbative calculation 
and in fact — > in the chiral limit, to^ 0. From 
the values above we get go = —82(11), g\ = —5.3(3) and 
52 = 77(12). On the other hand, the double scattering 
contribution to the ttq? s-wave [l8|, yields the identifica- 
tion gi — go — 92 — Ma{r^^)al P which has the correct 
pion mass dependence, provided all quantities are evalu- 
ated in the chiral limit, and suggests that 51—50 — 52 > 0. 



For 



0.5fm, a realistic value, and using either / 



or yields 51—50 — 52 ~ 40 — 47, while we get ~ — 20 
instead. A expected decrease of {r~^)a in the chiral hmit 
might naturally accomodate the discrepancy. On top of 
this, pion loop corrections to 7ra scattering, which are 
0(1//^), might have a sizable impact on 50,51 and 52 
due to chiral logs. Clearly, a more systematic assessment 
of the input values of the couplings and their uncertain- 
ties would be most useful. 

Let us now turn to the calculation of the long distance 
aa potential. The leading direct t-channel TPE contri- 
bution is depicted in Fig. [1] and can be written as 



K TVa — ^TTQ I 



{2i:Y{{p-qY~ml){p'^-ml) 



(4) 



There is a factor 5°''^6ab = 3 coming from tt'^tt" as well 
as 7r°7r° exchange (we neglect here tiny isospin breaking 
effects). As we see by power counting the integral is 
quartically divergent. Using the dispersion relation 
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with p(p^) the spectral density and c.t. stands for coun- 
terterms, e.g. A in Eq. ([T]), which will not contribute 
to the coordinate space potential at positive but non- 
vanishing distances. From Cutkosky's rules one gets 
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where A± = (51 — 5o)m^ -|-5o/i^/2±52(^^/4 — m^). Note 
that pf^^) is a positive quantity. This result agrees with 
Ref. [8| only when 52 = and resembles a similar calcu- 
lation for the central NN force [l3|. In the heavy mass 
limit ^/s — *■ 2Ma and in the CM system v ■ q = so that 
the potential is given by the expression 
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Computing the q and the /x integrals, up to distributions 
located at a; = 0, the final result becomes 
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iml [Ko{2x)f{x)+Kx{2x)g{x)] 
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(8) 



where x = m-^r, K{){x) and Ki{x) are modified Bessel 
functions and f{x) — 4(50 + 5i)'^a;** + 10(125q -|- 45250 + 
35i)+ [845o + 24(51 + 52)50 + 52(451 + I552)] x^, g{x) = 
4(5o+5i)(65o+52)a:3 + 10(1252-h45o52 + 35|)a:^. The TPE 
potential is attractive everywhere as can be recognized 
from the positivity of the spectral function. The potential 
^ is the strong interaction analog of the time-honoured 
Casimir-Polder electromagnetic forces [2^, [2l|. In fact, 
the TPE potential becomes singular at short distances, 



15(125^ + 45052 + 35!) 1 
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This is a relativistic and attractive Van der Waals inter- 
action which is exphcitly independent on the pion mass. 
In the opposite hmit of long distances we have 



(MeV) 



3(go+gi)^mre-^"^''' 
167r5/2M2/4 r5/2 



(10) 



Of course, the previous potential describes the strong 
interaction piece, and as we see is 0(1//*). Electromag- 
netic effects correspond to minimally couple photons. We 
keep one photon exchange (see Fig. [T]) but neglect two or 
higher photon exchanges as well as terms (^(Z^e^//^) 
which may be systematically computed from higher di- 
mensional corrections to the Lagrangean ([T]). These and 
other effects will be analyzed in more detail elsewhere. 

The total potential in a long distance expansion is 
given by adding the TPE and Coulomb potentials 



V{x) = V^l{x) + 



(11) 



with = 2 and = 1/137.04 the fine structure con- 
stant. The dots in Eq. ([TT|) represent shorter range cor- 
rections than TPE. Rotational invariance allows to write 
the relative s-wave function as ^(a?) = uo,p('')/\/47r with 
uo.p{r) the reduced s-wave function fulfilling 
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with U^lir) = Af„F„2-(r), o b = 2 /(M^Zge^) = 3.63fm 
the Bohr radius and p — \/MaE the CM momentum. 
The problem is to solve Eq. with suitable bound- 
ary conditions, but since the potential is singular and 
attractive at the origin, see Eq. ([S]), some renormaliza- 
tion proves necessary. Actually, the regularity condition, 
Wo,p(0) =0, only fixes the solution up to an arbitrary 
constant [23 . [23j which must be fixed independently of 
the potential. Furthermore, from the self-adjoint con - 
dition(for a discussion within the NN context see 
there is the relation 



21) 



= [uq^k 



(13) 



where the limit Tc ^ for the short distance cut-off is 
understood. Taking p — k and fc = we get further 
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(14) 



This means that the logarithmic derivative at the origin 
is energy independent and real leaving only one free pa- 
rameter left which will be fixed below. In essence, this is 
the non-perturbative renormalization program with one 
counterterm described in ^] for singular potentials. 

The *Be nucleus in its (J^,r) = (0+,0) ground state 
is unstable against a-decay, ®Be ^He -I- ^He with 
Q = 91.84 ± 0.04KeV and a very smaU (Breit-Wigner) 
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FIG. 2: The / = total two-pion exchange plus Coulomb po- 
tential barrier for the a — a system (in MeV) as a function of 
the relative distance (in fm). The horizontal straight line rep- 
resents the energy corresponding to the energy Q — 91.84KeV 
of *Be nucleus in its ground state, = 0'*'. 



width Tbw = 5.57± 0.25eV (see [2 5j for a review). The 
relative CM momentum is p = \/2iiaaQ — 19.2MeV 
{2pLaa = Ma). The Corresponding de Broglie wavelength, 
A ~ lOfm, is much larger than the size of the a particle, 
so one would not expect internal structure playing a cru- 
cial role. The outer classical turning point is determined 
by the Coulomb potential yielding rmax = 62.6fm. The 
situation is illustrated in Fig. [5] where the total I = 
TPE plus the Coulomb potential barrier for the a — a 
system are depicted, together with the experimental Q- 
value Q = 91.84KeV for the ^Bc "^Hc + "^He reaction 
which proceeds by standard tunnel effect. 

It is instructive to analyze the decay process within the 
WKB approximation where the width is given by ,26j 
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g-2 /,r».'^" dr\p{r)\ 



(15) 
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with p{r) = ^2[iaa[Q - V(r)). Here, 
are the classical turning points fulfilling V(riiiin) — 
l^(rmax) = Q = MsBe - 2Ma. Wc obtain FwKB = 
8.6(4)eV for the experimental Q, a correct order of mag- 
nitude compatible with the expected accuracy of the 
WKB formula. Although the TPE potential diverges at 
short distances, the inner classical turning point takes 
typically the value rrmin = 3fm for which there is about 
IMeV cancellation between TPE and Coulomb poten- 
tials. So, the tunneling region is not determined by 
the singularity. Finite cut-off corrections to Fwkb are 
0(rc^^) for Tc ^ Tinin as can be seen from Eq. (fT5| . 

A rigurous treatment of ^Be as a exponentially time- 
decaying state requires finding a pole of the S-matrix in 
the second Riemann sheet of the complex energy plane, 
so we look for exact numerical solutions of Eq. ful- 
filling the asymptotic boundary condition of a spherically 
outgoing Coulomb wave. 
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with ri = l/(pas) and p = pr. For complex momenta 
P = Pr + ipi the energy also becomes complex E = Q — 
iV 12. The boundary condition, Eq. ([H)) . implementing 
sclf-adjointncss provides a correlation between F and Q 
through the TPE potential. We get Tpoio = 3.4(2) eV 
for the S-matrix pole width, fairly independently of the 
cut-off radius for rc ^ r^i^ ~ 3fm. 

The experimentally determined [2^ Breit-Wigner 
small width involves the s-wave phase shift 27, [isj l 
Few = "i/SoiER) for 5o{Er) = Tr/2. We get 

rBw(^Be^aa) =4.3(3) eV, (exp.5.57(25)eV) , (17) 

for En — 91.8KeV and the couplings go, gi and 32 with 
their uncertainties obtained from low energy ira scatter- 
ing, Eq. ([3]) and Further, we analyze the scattering 
length ao and the effective range tq defined from the s- 
wave phase shift low energy expansion 

27rcot5o(p) 2 , , , 1 1 9 , „x 

f 2°. + —HV = + i,r,p^ + ■■■ 18 

aB(e^'"' - 1) as ao 2 

with h(x) the Landau-Smorodinsky function [27|. ^ 
From the universal low energy theorem of Ref. 
the Coulomb case we obtain (in fm) 
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1.03(1) 



5.3(3) 29(4) 



(19) 



The numerical coefficients depend on the total long dis- 
tance potential p4|) only. Using Eq. (fH)) we find a^^ — 
— 1210(70)fm and = 1.03(l)fm, in reasonable agree- 
ment with ao — — 1630(150)fm and ro — 1.08(l)fm from 
a low energy analysis of the data ^27]. The correla- 
tion between (Q,r) and (aoj^o) was pointed out long 
ago ^28,, 2^ but has no predictive power. As we see, 
this is compatible with the underlying chiral TPE poten- 
tial which, in addition, correlates tq with ao and F with 
Q. At higher energies further ingredients are needed; for 
Ei^AB = 3MeV, i.e. 1/p - 2.5fm, we get So = 141 ± 2° 
whereas Sq — 128.4±1° from the scattering data [2]. Note 
that we have made no attempt to refit the couplings go, 
gi and 32 stemming from low energy ira scattering, al- 
though it is obvious that further attraction is needed (a 
20 — 30% increase in go would do). The fact that we 
produce numbers in the right order of magnitude is en- 
couraging and calls for further improvements. 

A more systematic determination of the input param- 
eters would be highly desirable, possibly including pion 
loops in 7ra scattering. Moreover, an accurate description 
of aa scattering data in the higher energy elastic region 
below p-|-^Li and n-|-^B production threshold or a-decays 
of the excited ^Be states becomes sensitive to deeper in- 
teraction regions. This may require some peripheral nu- 
clear structure information involving e.g. a t + p and 
a He -|- n intermediate state contributions on top of 
the chiral interactions deduced here. With these provisos 



in mind, the present work provides an example on how 
chiral symmetry might be useful in low energy nuclear 
reactions eventually providing theoretical constraints in 
the hardly accessible regions of astrophysical interest. 
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